npj Breast Cancer

BackgroundBreast cancer exhibits extensive molecular heterogeneity across intrinsic subtypes, yet convergent metabolic reprogramming may represent an obligate feature of tumour initiation. We hypothesised that suppression of nuclear-encoded mitochondrial fatty acid oxidation (FAO) constitutes such a convergence point, defining a shared metabolic phenotype independent of subtype. MethodsRNA-seq data from 1,106 primary breast tumours and 113 normal-adjacent tissues (TCGA-BRCA) were intersected with 1,079 nuclear-encoded mitochondrial genes from MitoCarta 3.0. Differential expression was assessed using Welch t-test with Benjamini-Hochberg correction at all tumour stages, at Stage I specifically, and stratified across PAM50 subtypes. A 55-gene core FAO signature was derived by three-way intersection. Ten candidate genes were selected by pre-specified objective scoring, locked before any clinical testing. Gene set enrichment analysis (GSEA) was performed using MitoCarta 3.0 pathway annotations. Diagnostic performance, clinical associations, survival, and mutation independence were characterised. External validation used two independent GEO cohorts (GSE42568, n = 121; GSE109169, n = 50); prognostic validation used METABRIC (Molecular Taxonomy of Breast Cancer International Consortium; n = 1,980). DESeq2 was applied as methodological cross-validation. ResultsAmong 126 differentially expressed mitochondrial genes, fatty acid oxidation was the most significantly depleted pathway (normalised enrichment score -2.130; false discovery rate 0.001). The 55-gene core signature replicated in both external cohorts with 100% directional concordance (hypergeometric p < 10-{superscript 1}). All 10 candidate genes discriminated tumour from normal tissue (area under the curve 0.915-0.979) and demonstrated broad clinical associations. The composite FAO suppression score predicted overall survival in METABRIC (log-rank p = 7.82 x 10-) and MAOA achieved independent prognostic significance in multivariable Cox regression (hazard ratio 0.890; adjusted p = 0.009). DESeq2 cross-validation confirmed Spearman {rho} = 0.980 concordance. ConclusionsNuclear-encoded FAO suppression is a robust, pan-subtype feature of breast cancer detectable at Stage I and validated across independent platforms and cohorts. These 10 candidate genes constitute a consistent initiation-phase mitochondrial signature, implicating FAO suppression as a potential convergence point in breast cancer oncogenesis and motivating targeted functional investigation.

Purpose: Beyond estrogen receptor (ER) positivity, no genomic biomarker reliably identifies ER+ breast cancer patients who derive differential benefit from endocrine therapy (ET). We performed an unbiased genomic screen to discover genes predicting ET response and characterized the top candidate across clinical settings, treatment modalities, and an independent validation cohort. Experimental Design: We screened 240 genes in 1,197 metastatic ET-treated patients from the MSK-CHORD clinical genomics database using Cox proportional hazards regression with false discovery rate (FDR) correction. The top candidate, core-binding factor subunit beta (CBFB), was characterized across four cohorts defined by disease setting (metastatic/adjuvant) and treatment (ET/chemotherapy), with multivariable adjustment, gene-by-treatment interaction testing, left-truncation sensitivity analysis for guarantee-time bias, and external validation in METABRIC (N = 1,499 ER+). Results: CBFB mutations (prevalence, ~5%) were the only gene associated with improved time to progression (TTP). In metastatic ET patients, CBFB-mutated tumors (n = 80) demonstrated significantly longer TTP (hazard ratio [HR], 0.44; 95% CI, 0.29-0.67; P = .0002, FDR q = .010) with no chemotherapy benefit (HR, 1.16; P = .65). The gene-by-treatment interaction was significant (HR, 0.37; P = .009). Effects were robust to multivariable adjustment (HR, 0.46-0.50), independent of histology, and preserved under left-truncated Cox regression (HR, 0.38). In the adjuvant setting, CBFB mutations predicted improved recurrence-free survival (HR, 0.52; 95% CI, 0.31-0.85; P = .010), with no effect under chemotherapy. In METABRIC, CBFB mutations predicted improved ER+ overall survival (HR, 0.52; P = 9.3e-5). Conclusions: CBFB mutations identify ~5% of ER+ breast cancers with exceptional ET benefit. As CBFB is included on all major cancer gene panels, this biomarker requires no additional testing infrastructure for clinical implementation.

PurposeTo develop and validate a multimodal recurrence-risk model integrating histology, genomic testing, and clinical variables. MethodsWe developed AI-Path, a whole-slide image biomarker for recurrence prediction trained in CALGB 9344, and validated it in three independent cohorts: TAILORx, a multi-site Chicago cohort, and the MDX-BRCA cohort. We then integrated AI-Path with Oncotype DX Recurrence Score (RS), tumor size, and nodal status into a Cox model, PathClinRS, fit using 60% of cases from TAILORx, with the remaining 40% held out for validation. The primary end point was distant recurrence-free interval. Performance was assessed using Harrells concordance index (C-index) and Kaplan-Meier analyses. ResultsA total of 12,418 patients were included. In TAILORx, AI-Path outperformed RS for distant recurrence (C-index, 0.682 vs 0.647; P = .038), driven by superior prediction of late recurrence (0.656 vs 0.567; P < .001). In node-negative disease, PathClinRS outperformed RSClin in the TAILORx fitting (0.72 vs 0.70; P = .016) and validation sets (0.74 vs 0.70; P = .004). In node-positive disease, PathClinRS outperformed RSClinN+ in Chicago (0.94 vs 0.74; P < .001) and MDX-BRCA (0.71 vs 0.66; P = .004) cohorts. Compared with NATALEE eligibility, PathClinRS identified nearly twice as many high-risk node-negative patients while maintaining a comparable 10-year distant recurrence risk (16.7% vs 16.6% per NATALEE eligibility in TAILORx fitting; 21.0% vs 19.4% in TAILORx validation). PathClinRS identified 68% of intermediate risk premenopausal patients as low-risk with no evidence of chemotherapy benefit, compared to only 36% identified as low risk by standard clinicopathologic criteria. ConclusionDigital histopathology provides prognostic information complementary to genomic assays and has the potential to personalize therapy beyond existing clinicogenomic tools.

PurposeTriple-negative breast cancer (TNBC) exhibits substantial clinical heterogeneity, with some patients experiencing early recurrence and poor survival despite similar clinicopathologic features. We sought to determine whether quantitative measures of intratumoral collagen architecture and composition derived from standard histopathologic specimens can identify patients at risk of recurrence and adverse survival outcomes. Experimental DesignWe analyzed a retrospective cohort of 79 TNBC tumors assembled into a tissue microarray using a multimodal computational pathology framework integrating Massons Trichrome staining with COL1 and COL3 immunohistochemistry. Collagen architecture was quantified using fiber-based image analysis and unsupervised clustering, while collagen composition was assessed using a normalized COL3:COL1 ratio. Associations with recurrence-free interval (RFI) and overall survival (OS) were evaluated using Kaplan-Meier analysis, restricted mean survival time (RMST), and Cox proportional hazards modeling. ResultsUnsupervised analysis identified four distinct collagen architectural states, which were consolidated into low-risk and high-risk groups based on recurrence patterns. High-risk collagen architecture was associated with significantly worse long-term RFI (log-rank p=0.025; RMST difference 10.1 months). Independently, a higher COL3:COL1 ratio was associated with improved OS (log-rank p=0.042; RMST difference 9.4 months). Integration of architectural and compositional biomarkers further refined risk stratification, identifying a subgroup with high-risk architecture and low COL3:COL1 ratio that exhibited the poorest survival outcomes. Notably, collagen-based stratification identified patients with divergent outcomes not readily predicted from tumor stage alone. ConclusionsQuantitative assessment of intratumoral collagen architecture and composition provides clinically meaningful prognostic information in TNBC and enables stratification of recurrence and survival risk. These findings support extracellular matrix phenotyping as a practical and scalable computational pathology approach for refining risk assessment in TNBC. Translational RelevanceTriple-negative breast cancer (TNBC) remains clinically challenging due to heterogeneous outcomes that are not fully captured by standard clinicopathologic variables. In this study, we demonstrate that quantitative features of intratumoral collagen architecture and composition, derived from routine pathology specimens, provide clinically meaningful prognostic information. Collagen-based biomarkers, including distinct collagen architectural phenotypes and the COL3:COL1 ratio, identify patient subgroups with distinct recurrence and survival outcomes, particularly among individuals whose risk is not adequately predicted by conventional staging. Importantly, these features can be extracted from widely available histological stains and immunohistochemistry, supporting the potential integration into existing pathology workflows. These findings support the tumor microenvironment as an underutilized source of biomarkers and suggest that extracellular matrix-based phenotyping may improve risk stratification and inform clinical decision-making in TNBC.

Introduction: With recent approvals of multiple targeted therapies for triple-negative breast cancer (TNBC), including antibody-drug conjugates and immunotherapy in biomarker-selected populations, it is critical to define the temporal evolution of cell-surface target expression from early-stage to metastatic disease, the co-expression patterns across these markers, and optimal quantification methodologies. Here we report biomarker expression profiles measured by multi-omics and pathology-based platforms in patients with TNBC using a large cohort of matched longitudinal tumor samples. Methods: Patients who underwent neoadjuvant chemotherapy (NAC) for stage I-III TNBC or were diagnosed with any stage TNBC and developed metastatic recurrence were retrospectively identified from an institutional database and prospective research metastatic biopsy protocol. Tumor samples from diagnosis (DX), residual disease (RD) post-NAC (if applicable), and metastasis/recurrence (MR) were collected. Quantification of HER2, TROP2, and PD-L1 expression was performed by immunohistochemistry (IHC), whole-exome sequencing, transcriptome sequencing, and targeted mass spectrometry (MS). For HER2, TROP2, and stromal tumor-infiltrating lymphocytes (sTILs), both manual pathologist assessment and computational pathology quantification were obtained. HER2 status was categorized as HER2-0 or HER2-low by local (L-IHC) and central (C-IHC) review, TROP2 status was defined as low (H-score <100), medium (H-score 100-200) or high (H-score >200), and PD-L1 as low (tumor area positivity, TAP <5%) or high (TAP [&ge;]5%). Pathologist-assessed sTILs were classified as low (<10%), medium ([&ge;]10% and <40%) or high ([&ge;]40%). Biomarkers were compared between primary (DX/RD) and MR, and between pre- vs post-NAC (DX-RD) samples. Correlations between markers, quantification methods, inferred PAM50 subtype, and clinical variables of interest were evaluated. Results: A total of 359 samples from 110 patients with TNBC with data available from at least one platform were included in the analysis. HER2-low prevalence at DX, RD, and MR was: 51% (50/98), 40% (21/53), and 27% (16/60); TROP2 high/medium was 90% (47/52), 91% (42/46), and 88% (28/32); PD-L1-high was 51%, 50%, and 38% (9/24); and sTILs-high/medium was 88% (59/67), 80% (40/50), and 49% (17/35), respectively. While TROP2-high/medium vs low remained stable over time, HER2 IHC and sTILs significantly decreased from DX/RD to MR samples, both at the cohort-level (HER2, p=0.0081; sTILs, p=4.6x10e-5) and longitudinal patient-level (HER2, p=0.030; sTILs, p=0.0077), with a similar decreasing trend for PD-L1 that did not reach statistical significance. HER2 concordance (0 vs low) between L-IHC and C-IHC was 78% (91/116). ERBB2, TACSTD2 and CD274 mRNA expression were significantly correlated with IHC protein levels, though only TACSTD2 had limited overlap in distribution of gene expression between high/medium vs low groups. Strong correlation between protein membrane staining intensity from computational pathology, protein expression measured by MS, and pathologist-assed IHC was observed across all biomarkers tested by each method. In comparisons between biomarkers, pathologist-assessed PD-L1 IHC and sTILs were significantly correlated (p=0.0001); 94% (51/54) of PD-L1-high tumors were classified as sTILs high/medium. PAM50 subtype was not significantly correlated with time point or biomarker status, although there was a trend toward more HER2-enriched tumors in HER2-low (20%, 5/25) vs HER2-0 (6%, 3/52) (p=0.086). Across biomarkers and clinical variables, an association between age and sTILs was observed (p=0.038, FDR=0.42) due to a decrease in sTILs high/medium tumors with age, primarily driven by post-treatment (RD/MR) but not DX samples. Conclusions: Multi-platform and multi-omics profiling in this large unique cohort of longitudinal TNBC samples revealed distinct patterns of expression and dynamic changes of key biomarkers of interest for targeted therapies. Given variability with manual IHC scoring, improved methods for quantification of expression may help optimize treatment selection in an individualized manner.

BackgroundThe ketogenic diet is being explored as an adjuvant intervention in breast cancer because it lowers circulating glucose and elevates ketone bodies such as {beta}-hydroxybutyrate (BHB), but how individual ER+ breast cancer subtypes adapt to these conditions remains poorly characterized. We examined metabolic responses to BHB supplementation under glucose restriction in two ER+ breast cancer cell lines, asking whether metabolic adaptation patterns differ between models. MethodsMCF-7 and T47D cells were cultured under high glucose, glucose-restricted (5% of standard), or glucose-restricted with 10 mM BHB conditions and profiled by comprehensive two-dimensional gas chromatography-mass spectrometry (GCxGC-MS). Pairwise Welchs t-tests with Benjamini-Hochberg false discovery rate (FDR) correction were applied to identify treatment-responsive metabolites. Targeted assays quantified intracellular glycine, SHMT1 protein, and total branched-chain amino acid (BCAA) concentrations across a BHB dose range (2.5-15 mM). Patient tumor transcriptomic data from TCGA (n=1,084) and paired tumor-normal samples from GSE58135 (n=20) were analyzed for genes involved in one-carbon, ketone body, and BCAA metabolism. ResultsMCF-7 and T47D cells exhibited markedly divergent metabolic responses to BHB. In MCF-7 cells, BHB supplementation produced a broad pattern-level metabolic shift: 75% of detected metabolites trended upward when BHB was added to glucose-restricted cultures (C vs. B comparison), with 1,4-butanediol reaching nominal significance (FC=2.35, p=0.016) and a 4.1-fold trend increase in lactic acid (p=0.11), although no individual metabolite survived FDR correction. T47D cells showed essentially no metabolic response to BHB at the global level. Targeted assays detected an elevation in glycine at 5 mM BHB in both cell lines that did not follow a monotonic dose response and was not accompanied by changes in SHMT1 protein expression. Total BCAA levels were elevated by BHB in T47D cells but remained unchanged in MCF-7 cells. In paired patient samples, OXCT1 (log2FC = -1.41), SHMT1 (log2FC = -1.31), and ACAT1 (log2FC = -1.07) were significantly downregulated in ER+ tumors relative to matched normal tissue (adjusted p < 0.001 for all three). ConclusionsER+ breast cancer cell lines show heterogeneous metabolic responses to BHB supplementation under glucose restriction. The broad pattern of metabolite elevation in MCF-7 but not T47D cells suggests that capacity to utilize ketone bodies as metabolic substrate varies between ER+ models. The downregulation of OXCT1, ACAT1, and SHMT1 in ER+ tumors compared to normal tissue identifies these enzymes as candidate biomarkers that may help stratify which patients are likely to benefit from ketogenic interventions. Findings related to individual metabolites should be regarded as exploratory and require validation in larger, adequately powered cohorts.

Purpose: The prognostic role of tumor-infiltrating lymphocytes in luminal breast cancer remains uncertain, partly because density-based metrics do not capture spatial interactions between immune cell subsets. We developed a density-independent spatial metric quantifying macrophage-T cell proximity and assessed its prognostic value. Experimental Design: Using multiplex immunohistochemistry across three breast cancer cohorts (exploratory, n = 17; discovery, n = 687; validation, n = 305), we measured nearest-neighbor distances from T cells to M1-like and M2-like macrophages, benchmarked against a randomly subsampled total macrophage pool. We defined the Macrophage Spatial Polarity Index (MSPI) as the difference between M2-to-T cell and M1-to-T cell affinity scores, where higher values reflect an M2-dominated spatial phenotype. Cox regression was used to assess associations with distant disease-free survival (discovery) and overall survival (validation). Results: M2-like macrophages preferentially localized near T cells, independent of cell density. Higher MSPI was associated with shorter survival in luminal cancers (discovery: HR = 1.45, p < 0.001), with the strongest effect in young women with early-stage disease (HR = 2.16, p < 0.0001). MSPI remained independently prognostic after adjustment for stage, systemic treatment, and diagnosis period (HR = 2.31, 95% CI 1.73-3.09, p < 0.0001) and was non-significant in HER2-positive and triple-negative subtypes. Validation in an independent ER-positive cohort confirmed the finding (HR = 1.30, p = 0.004). Pooled analysis yielded HR = 2.13 (95% CI 1.68-2.70, p = 3.45 x 10-10). Conclusions: MSPI is a robust prognostic biomarker in luminal breast cancer, particularly in young women with early-stage disease, warranting further validation for risk stratification and therapeutic guidance.

BackgroundTriple-negative breast cancer (TNBC) exhibits substantial molecular heterogeneity and lacks targeted receptor therapies. Single-omic approaches inadequately capture its regulatory complexity, necessitating integrated multi-omic frameworks to identify stable prognostic signatures. MethodsMatched transcriptomic and DNA methylation data from the TCGA-BRCA cohort were normalised and mathematically integrated to isolate disease-associated variations. A calibrated machine learning voting ensemble (comprising LightGBM, Random Forest, and Logistic Regression) was trained to predict clinical survival. Model generalisability was tested on an independent microarray cohort (GSE58812) using independent quantile normalisation. SHAP (SHapley Additive exPlanations) values provided biological interpretability. ResultsDifferential and integrative analyses identified a 47-gene master prognostic signature. The ensemble classifier achieved an external validation accuracy of 74.77% (AUC 0.590) on unseen clinical patients. SHAP analysis confirmed the biological directionality of these specific biomarkers in driving mortality. Hypergeometric pathway enrichment highlighted targetable metabolic and signalling networks. ConclusionsThis multi-omic machine learning pipeline identifies a highly prognostic 47-gene signature for TNBC. The model demonstrates strong cross-platform generalisability and offers interpretable clinical utility for stratifying patient risk and guiding future therapeutic target development.

INTRODUCIONTumor cell infiltration in regional lymph nodes is a strong prognostic marker, guiding treatment decisions in breast cancer. While the immune cell composition in primary tumors has been more widely explored in later years, the immune cell composition of the sentinel node (SN) and axillary lymph nodes (ALN) remains understudied. A better understanding of how primary tumor and metastatic tumor cells alter the nodal immune microenvironment can shed light on metastasis and cancer progression to unveil new treatment strategies. MATERIALS AND METHODSFrom a prospective clinical cohort of 458 treatment-naive patients with primary operable breast cancer, we performed comprehensive immunophenotypic analysis using mass cytometry analysis of non-metastatic (SN-) and metastatic (SN+) and ALN (ALN+) lymph nodes. RESULTSAs expected, patients with ALN+ cases had a shorter time to distant metastases than SN+ and SN- cases. We identified an exhausted T-cell phenotype and an increase in Germinal Center B (GC B) cells and plasma cells in ALN+ samples compared to SN- samples, both in the whole cohort as well as when investigating estrogen-receptor positive (ER+) patients only. There were no differences in immune cell composition across breast cancer (BC) subtypes within SN-samples. SN+ samples from triple negative BC (TNBC) showed a trend towards increased abundance of GC B and plasma cells, similar to more advanced ALN+, suggesting that smaller TN metastases may trigger an immune activation at an early stage of dissemination. Further analysis of SN- samples from ER+ patients revealed a subset of patients where the immune response had a more exhausted T-cell phenotype. This group was enriched for lymph nodes that were deemed negative by ordinary pathology examination (microscopy) but had detectable tumor cells by CyTOF analysis. CONCLUSIONThe immune profiles of SN and ALN samples from breast cancer patients are highly diverse, showing limited associations to BC subtype, clinical parameters or patient outcome. Metastatic tumor cells play a significant role in driving T-cell exhaustion and immunosuppression. Notably, in approximately 50% of the ER+ samples, T-cell exhaustion was detectable. This coincides with the presence of tumor cells identified by CyTOF, which were likely missed by conventional pathological examination. These findings suggest that small tumor deposits alter the immune composition, and the immune profile reveals the presence of tumor cells.

BackgroundRecent evidence suggests that cancer can exhibit splicing alterations that give rise to tumour-specific isoforms. One example is NUMB, which produces four isoforms (p72, p71, p66, and p65) through alternative splicing of exons 3 and 9. Traditionally considered a tumour suppressor, it also has been considered an oncogene. We propose that this duality is due to isoform-specific expression. ResultsUsing public databases, we identified a tumour-associated switch in NUMB isoform expression: p72/p71 are upregulated in tumours, whereas p66/p65 are more expressed in non-tumour tissues. These isoforms correlate differently with cellular processes. NUMBL, a NUMB homolog, behaves similarly to p65. We identified two transcriptional clusters: one characterized by high expression of p72/p71, and another by p66/p65/NUMBL. Each group was associated differently with the Notch, WNT/{beta}-catenin, Hedgehog, and Hippo signalling pathways, suggesting isoform-specific regulatory roles. Analysis of breast cancer cell lines (CCLE) led to a NUMB score based on isoform expression, which classified cell lines into biologically distinct groups. The p72/p71-enriched group showed distinct signatures, pathway activity, and drug sensitivity. Applying this score to TCGA-BRCA samples revealed a significant link between high NUMB-score and poor survival, confirmed by Kaplan-Meier analysis. ConclusionsNUMB emerges as a potential oncogenic contributor and biomarker in splicing-based personalised medicine, highlighting isoform-specific expression as a clinically relevant determinant of tumour behaviour, pathway activity, and therapeutic response.

Assessment of pathologic complete response (pCR) following neoadjuvant chemotherapy (NAC) remains an unmet clinical need in breast cancer. Fibroblast activation protein inhibitor (FAPI) PET targets the tumor microenvironment and may therefore enhance response evaluation after NAC. This study aimed to compare the performance of [68Ga]Ga-FAPI-04 PET, [18F]FDG PET, and contrast-enhanced MRI for predicting pathologic response after NAC in breast cancer, with separate analyses for primary breast lesions and axillary lymph nodes. MethodsIn this prospective single-center diagnostic accuracy study, women with biopsy-confirmed stage II-III breast cancer underwent baseline and post-therapy [68Ga]Ga-FAPI-04 PET/MRI, [18F]FDG PET/CT, and contrast-enhanced MRI before surgery. Quantitative PET parameters were evaluated for primary tumors and axillary lymph nodes. pCR was defined as ypT0/isN0. Significant variables identified in univariable analyses were further explored using least absolute shrinkage and selection operator (LASSO) analysis, and receiver-operating-characteristic (ROC) analysis was performed to assess diagnostic performance. Fibroblast activation protein expression was also assessed by immunohistochemistry in paired pre- and post-therapy tumor specimens from a subset of patients. ResultsTwenty-four patients completed the study protocol, yielding 25 primary lesions and 44 metastatic lymph nodes across 27 axillary compartments. Overall patient-level pCR was achieved in 13 of 24 patients (54.17%). The lesion-level pCR rate was 60.00% (15/25) for primary breast lesions, and the node-level pCR rate was 72.73% (32/44) for axillary lymph nodes. For primary tumor response, post-therapy [68Ga]Ga-FAPI-04 SUVmax showed the highest diagnostic performance (AUC, 0.84; sensitivity, 80.00%; specificity, 80.00%; accuracy, 80.00%), whereas the optimal [18F]FDG parameter was {Delta} TBR% (AUC, 0.747). For nodal response, post-therapy [68Ga]Ga-FAPI-04 SULmean showed the highest diagnostic performance (AUC, 0.89; sensitivity, 91.67%; specificity, 81.25%; accuracy, 84.09%) and was significantly different from the best [18F]FDG parameter ({Delta} SULmax%, AUC, 0.669) on DeLong testing (P < 0.05). MRI achieved AUCs of 0.733 for primary lesions and 0.770 for lymph nodes. Stromal FAP expression positively correlated with [68Ga]Ga-FAPI-04 SUVmax and was markedly reduced in lesions achieving pCR. ConclusionPost-therapy [68Ga]Ga-FAPI-04 PET may serve as a promising adjunctive imaging biomarker for predicting pathologic response after NAC in breast cancer, particularly for axillary nodal assessment. These findings suggest that FAPI PET may provide clinically relevant information for preoperative evaluation of residual disease burden, potentially contributing to more individualized surgical planning and treatment decision-making.

Purpose: Unilateral breast cancer (BC) patients scheduled for mastectomy often choose to undergo contralateral prophylactic mastectomy (CPM), despite substantial declines in contralateral breast cancer (CBC) risk in recent decades. Models predicting absolute risk of future CBC can aid informed decision-making about CPM. CBCRisk is an existing CBC absolute risk prediction model trained on unilateral BC patients regardless of whether they had mastectomy. Here we developed CBCRisk-Mastectomy, tailored specifically to BC patients scheduled for mastectomy and considering CPM. Patients and Methods: We used data on BC patients who underwent mastectomy to treat their first BC from two nationally representative sources: Breast Cancer Surveillance Consortium (BCSC) and Surveillance, Epidemiology, and End Results (SEER) cancer registry. We imputed missing data in the BCSC sample and used conditional logistic regression models, trained on 2,660 BC patients (665 CBC cases) from BCSC, to identify predictors and estimate relative risks (RRs). These were combined with attributable risks and CBC incidence rates estimated from SEER to obtain absolute risk. Cross-validation was used to internally validate CBCRisk-Mastectomy and compare with CBCRisk. Results: CBCRisk-Mastectomy has nine predictors: first BC type, lobular carcinoma in situ status, estrogen receptor status, tumor stage, breast density, age at BC diagnosis, family history of BC, age at first birth, and body mass index. The areas under the curve and their 95% confidence intervals for 5-year predictions for CBCRisk-Mastectomy and CBCRisk were 0.62 (0.59, 0.65) and 0.58 (0.55, 0.61), respectively. Conclusions: CBCRisk-Mastectomy may aid clinicians in counseling BC patients scheduled for mastectomy, enabling improved decision-making regarding CPM.

Chemotherapeutic treatment of breast cancer with Doxorubicin (DOX) can induce tumor and stromal cell senescence leading to therapy-resistance. Senescence-associated secretory phenotype (SASP) promotes secretion of pro-inflammatory and tumorigenic factors causing systemic inflammation. Combined, this can result in immune suppression, tumor growth and secondary spread of cancer. Targeting and removing senescent and cancerous cells using a combination of chemotherapeutic and senolytic drugs may reduce systemic inflammation, improve therapeutic efficacy, and prevent metastasis. Exposure of triple-negative breast cancer (MDA-MB-231), hormone-responsive (MCF-7) and HER2+ (MDA-MB-453) cells, and primary spine osteoblasts to DOX showed significant induction of p21-positive senescent cells. DOX and senolytics (RG-7112, o-Vanillin) treatment of co-culture spheroids showed a significant additive effect in reducing tumor sphere viability and growth, indicating reduced metastatic potential. This was correlated with reduced SASP in triple-negative and hormone responsive lines and decreased levels of senescent cells in all subtypes and primary stromal cells, while proliferation was decreased, and apoptosis increased across all breast cancer subtypes. Future chemotherapeutic treatment in breast cancer models may be optimized by adding senolytic drugs to more effectively clear senescent tumor and stromal cells, reducing risk for relapse and metastatic potential, while allowing for tissue regeneration in the bone metastatic environment. Graphical Abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=113 SRC="FIGDIR/small/724653v1_ufig1.gif" ALT="Figure 1"> View larger version (24K): org.highwire.dtl.DTLVardef@c4cb8forg.highwire.dtl.DTLVardef@105219org.highwire.dtl.DTLVardef@17e0517org.highwire.dtl.DTLVardef@802bd2_HPS_FORMAT_FIGEXP M_FIG C_FIG Senolytics selectively eliminate senescent cancer and stromal cells and enhance Doxorubicin efficacy in a 3D bone-like tumor microenvironment model.

Background Background breast features are frequently noted in pathology reports alongside invasive breast cancer but rarely factor into prognosis or treatment decisions. Their relationship to tumor characteristics and patient outcomes remains incompletely characterised. Methods We conducted a retrospective cohort study of 7,603 patients with Stage I-III invasive breast cancer (diagnosed 1991-2022, age <80 years) from the Joint Breast Cancer Registry in Singapore. Natural language processing (NLP) was applied to 9,754 free-text pathology reports to extract co-existing background breast features, with accuracy validated by dual-reviewer assessment of 200 reports. Unsupervised hierarchical clustering grouped extracted features into three categories. Associations with tumor characteristics were assessed by multinomial logistic regression, and ten-year overall survival by Cox proportional hazards models (median follow-up 9.6 years; 620 deaths). Results Here we show that NLP-based extraction of background breast features from routine pathology reports achieves an accuracy of over 90% across features. Lobular neoplasia and benign proliferative changes are associated with less aggressive tumor characteristics, whereas early neoplastic and papillary lesions are more prevalent in HER2-enriched and luminal B tumor subtypes. Benign proliferative changes are associated with better survival in age- and year-adjusted models (hazard ratio 0.91, 95% CI 0.86-0.97), but this association is attenuated after adjustment for stage and subtype. Conclusions NLP-enabled extraction of background breast features from pathology text is feasible at scale. These features reflect tumor biology but do not independently add prognostic information beyond established clinical variables.

Invasive lobular carcinoma (ILC) is the most frequently diagnosed special histological subtype of invasive breast cancer and accounts for 10 - 15% of all cases. The pathognomonic hallmark of ILC is the genetic loss of E-cadherin (CDH1) causing the disruption of adherens junctions and resulting in discohesive, linear growth. To better understand the role of E-cadherin in ILC metastasis, we generated three ILC cell lines, MDA-MB-134-VI, SUM44PE, and BCK4, with inducible E-cadherin expression, resulting in successful restoration of functional adherens junctions. E-cadherin expression reduced growth in 2D culture, and that effect was even greater in 3D ultra-low attachment (ULA) conditions where increased cell death was consistent with the previously described role of E-cadherin in anoikis. E-cadherin expression did not rescue the lack of migration and invasion of ILC cell line models; however, it decreased haptotaxis and increased adherence to Collagen I in SUM44 cells. There was no significant effect of E-cadherin expression on primary orthotopic tumor growth, but spontaneous metastasis to the reproductive tract, brain, and GI tract was reduced. Inhibition of metastasis to the reproductive tract and brain was also seen after tail vein injection of MDA-MB-134 E-cadherin-expressing cells. In summary, overexpression of functional E-cadherin in ILC models has some, but limited, effects on 2D growth in vitro and primary tumor growth in vivo, but there are pronounced effects on 3D ULA growth and metastases in vivo, with stronger effects on metastatic sites enriched in patients with ILC, especially the reproductive and GI tracts.

Background: Population-relevant BRCA1/BRCA2 data from Bangladesh are scarce, creating challenges for hereditary breast and ovarian cancer variant interpretation, counseling, and follow-up testing. We examined a clinically referred Bangladeshi cohort to characterize assay-derived BRCA1/BRCA2 short variants, sequencing-depth performance, and copy-number findings in a conservative pilot framework. Methods: Twenty-three de-identified blood-derived DNA samples were assessed using a targeted BRCA1/BRCA2 next-generation sequencing workflow. Downstream analysis used assay-generated short-variant, coverage, and CNV outputs, with coordinates reported on hg19/GRCh37. Short variants were evaluated from high-confidence PASS/VCC-H calls, and CNV review incorporated both target-region and amplicon-level copy-number patterns. Results: After removal of four low-VAF review observations, the primary germline-compatible dataset comprised 304 short-variant observations representing 34 unique variants. Both BRCA1 and BRCA2 contributed comparable variant burdens, while the overall profile was mainly composed of missense and synonymous changes. Six sample-specific heterozygous BRCA1 truncating candidates were observed, including five frameshift variants and one stop-gain variant. Protein-level mapping placed these events across the central-to-C-terminal portion of BRCA1. Sequencing depth was consistently high across the targeted regions, with all 4,255 amplicon-sample measurements exceeding 280x and 99.91% reaching at least 500x. Copy-number analysis highlighted one candidate BRCA1 multi-exon deletion-like event involving exons 15-20 in BCSIR-BRCA-21, with unresolved partial exon 14 involvement. Conclusions: This study provides an initial Bangladesh-focused targeted BRCA1/BRCA2 dataset and identifies candidate short-variant and CNV findings for validation. These findings should be interpreted as analytical candidates only and require confirmatory testing and expert clinical curation before any clinical application. The cohort is referral-enriched and should not be used to infer population prevalence.

Cancer-related fatigue is the most common and persistent symptom in breast cancer, with fatigue reported up to 10 years post-diagnosis. Unlike many cancers, fatigue in breast cancer often arises during early-stage disease in the absence of cachexia. While many factors contribute to fatigue, the direct contribution of cancer-associated skeletal muscle pathology remains poorly understood. Here we analyzed pectoralis major muscle biopsies from individuals with non-metastatic breast cancer and non-cancer controls. Using single-fiber morphometry and spatial transcriptomics, we identified fiber-type-specific structural alterations and spatially localized transcriptional reprogramming within the muscle microenvironment. Single-fiber morphometry revealed selective atrophy of oxidative type I and type IIa muscle fibers, while glycolytic type IIx fibers were relatively preserved. Concordant spatial transcriptomic profiling revealed suppression of oxidative metabolic programs, evidence of mitochondrial dysfunction, and spatially localized catabolic signaling originating from intramuscular adipocytes. This study introduces an integrated framework for profiling skeletal muscle architecture and spatially localized gene expression in surgically obtained muscle biopsies and represents the first application of spatial transcriptomics to human skeletal muscle from individuals with cancer. These findings demonstrate structural and metabolic remodeling of skeletal muscle in non-metastatic breast cancer and suggest targeting muscle metabolism represents a promising therapeutic strategy for cancer-related fatigue.

Purpose: Chemotherapy-induced peripheral neuropathy (CIPN) affects many people receiving taxane treatment for breast cancer. Symptom trajectories vary, with some recovering, and others experiencing persistent, or delayed worsening (coasting) symptoms. The prevalence and predictors of these trajectories remain unclear. This study identified the prevalence and biopsychosocial predictors of CIPN persistence, improvement, and coasting within three months post-treatment. Methods: This secondary analysis included participants treated with taxanes for stage I-III breast cancer who completed the Functional Assessment of Cancer Therapy/Gynecologic Oncology Group-Neurotoxicity-4 (FACT/GOG-NTX-4) at baseline, post-chemotherapy, and three months later. A minimally important difference (MID) from baseline on the FACT/GOG-NTX-4 defined persistence, improvement, coasting, and no MID-CIPN (below the MID threshold at each assessment) trajectories. Baseline assessments included self-reported pain/well-being, sensory, balance, and lower limb physical functioning measures, and sociodemographic and treatment data were collected. Results: Among 102 participants (51.57{+/-}11.24 years), persistence occurred in 34.3%, improvement in 25.5%, coasting in 6.9%, and no MID-CIPN in 33.3%. Compared to no MID-CIPN, older age (OR=1.120; 95%CI: 1.026-1.222), higher expected pain (OR=1.630; 95%CI: 1.082-2.456), and cold hyperalgesia at the foot (OR=1.130; 95%CI: 1.018-1.254) predicted persistence. Lower fatigue predicted improvement (OR=0.904; 95%CI: 0.845-0.968). No predictors were identified for coasting. Conclusion: CIPN trajectories are heterogeneous. Age and pre-treatment pain expectations, cold hyperalgesia, and fatigue differentiate patients with persistent CIPN and those likely to improve from those with no CIPN. Implications for Cancer Survivors: Early identification of individuals at risk for persistent neurotoxicity may support risk stratification and guide targeted supportive care strategies.

BackgroundThe MDM2-p53 signaling pathway plays a central role in tumor suppression, and genetic variants that disrupt this pathway may influence breast cancer (BC) susceptibility. However, data from South Asian populations, particularly Bangladesh, remain limited. MethodsA case-control study was conducted in Bangladeshi women, including BC patients and healthy controls (HCs). Genotyping of MDM2 polymorphisms was performed using PCR-based methods. Circulating MDM2 and p53 protein levels were measured using enzyme-linked immunosorbent assays (ELISA). Associations between genotype, protein levels, BC status, and clinicopathological features were evaluated using appropriate statistical models. ResultsA strong and genotype-specific association was observed for MDM2 rs2279744. Women carrying the heterozygous TG genotype had a markedly increased risk of BC across additive, dominant, and over-dominant models, whereas the GG genotype showed a protective effect under the recessive model. In contrast, rs937282 did not show a significant association with BC risk. Circulating MDM2 levels were significantly elevated in patients compared with controls and varied by rs2279744 genotype, while circulating p53 levels showed an opposite trend. A strong inverse correlation was observed between serum MDM2 and p53 levels, supporting dysregulation of the MDM2-p53 feedback loop. Elevated MDM2 levels were also noted in HER2-positive and triple-positive BC subtypes. ConclusionTogether, these findings indicate that the MDM2 rs2279744 polymorphism contributes to BC susceptibility in a genotype-specific manner, likely through disruption of the MDM2-p53 regulatory balance. However, the absence of functional validation limits direct causal inference.

Quantitative biophysical signatures of nuclear spatial reorganisation across breast carcinoma progression remain insufficiently characterised. We apply two complementary fractal descriptors, Correlation dimension (Dc) and Minkowski dimension (Dm), to 4276 regions of interest across seven breast tissue subtypes from the BRACS dataset, validating observed dimensions against systematically constructed null spatial models to distinguish genuine structural organisation from geometric irregularity. All subtypes significantly exceed the complete spatial randomness baseline, confirming universal departure from random nuclear arrangement. The observed scaling is characterised as statistically monofractal within a bounded pre-fractal range. Invasive carcinoma uniquely fails to exceed the clustered null in Dc while simultaneously showing the weakest Dm null deviation, a dual convergence toward stochastic baselines consistent with the progressive removal of architectural constraints. Flat epithelial atypia exhibits a unique directional dissociation with the lowest Dc across all subtypes combined with high Dm null deviation, a co-occurrence not observed in any other subtype and geometrically consistent with decoupled nuclear spatial organisation at the centroid distribution and boundary morphology scales. Interpreted within a percolation-theoretic framework, the non-monotonic null deviation trajectory maps onto qualitative regime transitions, providing a physically grounded explanation for the observed discrimination profile across pathological transitions. These findings position fractal-like nuclear architecture as a potential descriptor for pre-malignant transitional states.